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The study on the charmed scalar meson spectroscopy has become a hot topic both experimen- 
tally and theoretically. The Bf s ) decays provide an ideal place to study their property. We employ 
the B-meson light-cone sum rules to compute the — > D* + (2317) and B~ — > Do" (2400) transi- 
tion form factors at large recoil, assuming D*+(2317) and Z?5 (2400) being scalar quark-anti -quark 
states. The results are extrapolated to the whole momentum region with the help of HQET. Con- 
sidering large uncertainties, our results can be consistent with the previous studies, while the power 
corrections should be large. We also estimate the semi-leptonic decays B® — > D* + (2317)^; and 
B~ -> L>o°(2400)ZPi. The branching fraction of the semi-leptonic B° ->• D*+(2317)Zz?i decay is 
around 6 x 10 -3 for light leptons and 0.8 x 10 -3 for tau final state. The predicted branching ration 
of B~ — s> DQ°(2400)i^i is slightly larger than B® — > D* + (2317)i^(, and we hope the future data in 
LHCb can test these results. 

PACS numbers: 14.40. Lb, 13.20.He, 11.55.Hx 



I. INTRODUCTION 

The charmed scalar meson spectroscopy has evoked many interests since the observation of _D* (2317) by Babar 
collaboration at 2003. In addition, the signal for the isospin doublet Dp (2400) has also been reported by Belle0and 
FocusQin the Dtt final state. Recently more measurements on the charmed scalar meson final state in B decays have 
been performed^. The low mass and the narrow width of -D* (2317) indicates some hints on its mysterious inner 
structure. It is regarded as a scalar meson state in some studies, while it has also been assigned to be a four-quark 
state or the molecular state. Until now, the structure of D* (2317) is still a controversial problem. As for Dq(2400), 
there is less information from experiments, and our knowledge of its property is even poorer. So we need more 
phcnomenological analysis to clarify the inner structure of these p-wave states. 

A great number of B decay events have been accumulated at B factories which provide good places to test the inner 
structure of the charmed scalar meson. To study the £?-to-scalar meson decay modes theoretically, an essential task is 
to evaluate B° -> L>+(2317) and B~ -> L>S°(2400) transition form factors. In heavy quark effective theory(HQET)[|, 
the heavy-to-heavy form factor can be reduced to the universal Isgur-Wise(IW) function £(i> • v') in the heavy quark 
limit. In order to estimate the form factors or the IW function, one must employ the non-perturbative methods. 
There have existed some phenomenological studies using different approaches, including the phenomenological model 
Q, the QCD sum rules approach 0-d, PQCD approach [l(|, Lattice QCD [IlT - [l3| . as well as the light-cone sum rules 

(LCSR) a 

LCSR [15H17] combines the traditional QCD sum rules [18[ with the theory of hard exclusive process, and offers 
a systematic way to compute the soft contribution to the transition form factor. The vacuum-to-hadron correlation 
function is computed in terms of light-cone OPE in the LCSR. In the conventional LCSR for B® — > Z)+(2317) form 
factor, the correlation function is taken between the vacuum and D+(2317) state, whereas the B meson is interpolated 
by a local current. The long distance effect of the form factor is then described by the distribution amplitudes(DAs) 
of D+(2317). As the structure of D+(2317) is not well understood, the DAs of £>+(2317) are rather model dependent. 
In this paper, we employ a different sum rule for the transition form factor following Ref.[l9|, where the correlation 
function is constructed with the on-shell .B-meson and the interpolated current for the charmed scalar meson. As the 
nonperturative dynamics is parameterized in terms of the B-mcson DAspol [2l| . the new method is usually called 
B-meson LCSR and it has been widely applied to the calculation of heavy-to-light matrix elements [22l. [23j . 

In this work, we will employ the .B-meson LCSR approach to evaluate the the B® -> D*^ (2317) and B~ -> 
Dq°(2400) form factors. In our calculation D*+(2317) and £>5°(2400) are regarded as qq mesonic states. The relevant 
semi-leptonic B® — > D*q(2317)Iu and B~ — > Dq°(2A00)Iv decay modes are also analyzed. The large number of data 
accumulated in the B factories and LHC-b can test whether our assumption is reasonable, and the result can help to 
clarify the inner structures of the new measured charmed scalar mesons. 

The paper is arranged as follows: We firstly derive the LCSR for the B" -> L>*+(2317) and B~ -> L>S (2400) 
form factors in the section II. The contributions from both two-particle and three-particle wave functions of B meson 
are computed. The numerical analysis of LCSR for the transition form factors at large recoil region is displayed in 
section III. The HQET is adopted to describe transitions at the small recoil region. Moreover, detailed comparisons 
between the form factors obtained under various approaches are also presented here. Utilizing these form factors, the 
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branching fractions of scmilcptonic decays are calculated in section IV. The last section is devoted to the conclusion. 



II. THE LIGHT-CONE SUM RULES FOR FORM FACTORS 



The B-to-charmcd scalar meson transition form factor induced by an axial vector current is defined by: 

(D* (p)\c W b\B(j> + q)) = -i{v»fU : M 2 ) + ^BD*(^} ' W 

where the notation denotes B°,B+ and B s , and Dq refers to D*^ (2317) and £>JS°(2400). To obtain the form 
factors with B meson LCSR, we consider the following correlation function with on-shell _B-meson state: 

F M (p,q)=i J d 4 xe tp - x (0\T{q(x)c(x),c(0)l^-l 5 )b(0)}\B(P + q)), (2) 

where C7 M (1 — 75)6 is the b — > c (electro)weak currents and qc is the interpolating current for a charmed scalar meson. 
The hadronic representation of the correlation function can be written as 

p / x (0|g(0)c(0) I J>g (p)) (gg (p) |c(0) 7 „7 5 b(0) \B(P + q)) 
F » M = rnl.-p> 

h y 

The decay constants fu* and fu* are given by 

(0\q^c\D*( P )) = f D . Plt , (0\qc\D*(p)) = mn S fn* , (4) 

where fo* = {m c — m q )f]j»/m]j* and m c , m q are the current masses of charm quark and light quark, respectively 
Inserting the definitions of the form factors and decay constants, the correlation function reads: 

— iwn* fn* 

p h + {s,q 2 )p IJ , + p h _{s,q 2 )q u 

as 5 ' ( 5 ) 

s — p z 

where Sq° is the threshold parameter corresponding to the Dq channel. 

On the other side, in the deep Euclidean region, the correlation function can be calculated in the perturbative 
theory using the operator production expansion near the light cone : 

F,(p,q) = F? CD ( q 2 ,p 2 )p, + F^(q 2 ,p 2 )q, (6) 

l ImFg CD (q 2 ,p 2 ) r llm^qV) 
as 5 pa + / ds t, g M . 

7r s - p z J m 2 7r s-p z 

Applying the quark-hadron duality 

P 1(s,q 2 ) = -lmF t 9 CD (q 2 ,p 2 )Q(s - s h ), (7) 

with i = "+, — " and performing Borel transformation with respect to the variable p 2 , we can derive the sum rules for 
the form factors as 
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The leading-order contribution to the OPE is illustrated in Fig. la. The correlation function can be calculated by 
contracting the charm quark fields in Eq. ^ and inserting the c quark propagator, then we arrive at: 

Fl B \p) = l Jd i x e*'J ^ie- ik - x (0\T{q(x)S F (x,0) 7v (l- 75 )b(0)}\B(P B )) (9) 
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FIG. 1: Diagrams corresponding to the contributions of (a) two-particle and (b) three-particle B -meson DA's to the correlation 
function flj) 



The full quark propagator can be written as[24|, 
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where the first term is the free-quark propagator and Gfj = G a pv T^ with Ti[T a T b ] = \8 ab . Inserting this propagator 
to Eq.®, we can find that the long distance contribution to the correlation function is expressed by non-local matrix 
elements, which defines the S-meson light-cone DA. In the leading Fock state: 
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pa 



where [a;, 0] is the path-ordered gauge factor. The variable uj > is the plus component of the spectator-quark 
momentum in the B meson. The three-particle DAs' contribution is shown in the diagram Fig. (lb), with the 
definition 
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(12) 



where the gauge link factors arc omitted for brevity. The DA's ^v^A, Xa and Ya depend on two variables w and 
£, corresponding to the plus components of the light-quark and gluon momenta in the B meson. 

Substituting the B meson distribution function into the correlation function and employing the quark hadron 
duality ([7]), we arrive at the sum rules for transition form factors as 
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where the argument of the wave functions is m^cr. In addition, a = 1 — a and <tq is the root of the equation 
as — (era + r 2 )m 2 B + om 2 B q 2 = 0. The modified wave function $±(w) = Jq cZt[</> + (t) — </>_(t)]. The contributions 

from three particle -B meson DAs are denoted by f BD ^' and /^zf* > which are given in the appendix. 



III. NUMERICAL ANALYSIS OF SUM RULES FOR FORM FACTORS 

Now we are going to calculate the form factors fjj* (q 2 ) and fjj* (q 2 ) numerically. In the following, we list the relevant 
input parameters for the D+(2317) and 1)^(2400). Their mass is taken from PDG [H]: m D ' = 2.318GeV and m D * = 
2.318GcV. The decay constant f D * o = (250 ± 25)MeV"[2fJ. For the £^(2400) state, we expect fb a /fb s0 = fb/ fb, in 
the SU(3) limit. We adopt the values f D = (223 ± 18)MeV and f Ds = (274±20)MeV, we find f D * = (203 ± 30)MeV. 
As for the decay constant of B s meson, we use the results fs = 130MeV [27| and JbJSb = 1.16 ±0.09 [28| determined 
from QCDSR. The threshold parameter sq can be fixed by fitting the LCSR of the charmed meson masses to the 
experimental data. Numerically, the threshold value in the X channel would be s x = {rax + Ax) 2 , where Ax is 
about 0.6 GeV j29Tj3l1 |. and we simply take it as (0.6 ± 0.1) GeV in the error analysis. The two-particle DAs of 
-B-meson inspired from QCD sum rule analysis reads [20j : 

Br \ u - — 
0+(w) = -2 e 
ui. 





B H = -e"^, (15) 

CJ 



and the 3-particle DAs are given by: 



6w5 

6wg 

Y A (u,t) - --^ I ^o-13c + 30e-^ + « ) /"°. (16) 
24wq 

The parameters ojo, \h and Xe satisfy the conditions adopted in(20l]: 

«* = §A, A| = A 2 h = ^ 2 = |a 2 . (17) 

Numerically we employ the values u)q = 0.45 ± O.lOGeV and ljq s = 0.50 ± O.lOGeV, here we have taken small SU(3) 
breaking effect into account. 

After fixing the corresponding parameters, we can proceed to compute the numerical values of the form factors. In 
principle, the form factors should not depend on the the unphysical Borel mass M 2 . However, the OPE series are 
truncated up to next to leading Fock state of the B meson and the QCD corrections are not considered, a manifest 
dependence of the form factors on the Borcl parameter M 2 would emerge. Therefore, we should search for the so-called 
"Borel window" , where Borel mass dependence is mild, in order that the truncation is acceptable. 

We firstly focus on the form factors at zero momentum transfer. To extract the form factor fb* (0), the contribution 
from the higher resonances and continuum states should be less than 30 % in the total sum rules and the value of 
fb* (0) should not be sensitive to the Borel mass. In view of these considerations, the Borel parameter M 2 should not 
be either too large or too small. To make sure that the contributions from the higher states are exponentially damped 
( see Eq. (fT4|)) and the global quark-hadron duality is satisfied, we need a smaller Borel mass. On the other hand, the 
Borel mass could not be too small for the validity of OPE near the light-cone for the correlation function, since the 
contributions of higher twist distribution amplitudes amount to the higher power of 1/M 2 to the perturbative part. 
In this way, we find a Borel platform M 2 G [3.5,5]GeV 2 . The Borel mass dependence of the form factors is plotted 
in Fig. [2] and Fig.[3j the former includes the contribution from the three-point B meson distribution amplitudes and 
the higher states contribution is shown in the latter one. From these diagrams we can easily see that the higher Fock 
state is highly suppressed in the Borel window, and higher exited states and the continuum states contribution is 
within 15% for f^, (0)(30% for /^»(0)). The numerical value for these form factors arc collected in Table |H where 
the uncertainties are from the combination of the variation of shape parameter loq, the fluctuation of threshold value, 
the uncertainties of quark masses and the errors of decay constants for the involved mesons. The results in the other 
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FIG. 2: The dependence form factor f£, (0) and f D » (0) on the Borel mass M 2 , the contribution from 2-point B meson DA 
is denoted by the dashed line, and the dotted line represents the 3-point DA contribution. The solid line gives the total results. 
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FIG. 3: The dependence form factor /^« (0) and /^« (0) on the Borel mass M 2 , the contribution of higher exited states and 
the continuum states in the whole sum rules is shown by the dotted line. 



studies are listed for comparison. We can see that for /^(O) our result is sightly larger than light-meson LQSR, 
however the results are consistent with each other within the errors. For f^, (0), the sign of our result is consistent 
with that obtained from the QCDSR, but it is different from that derived in the light meson LCSR. This discrepancy 
is expected to be smeared by power corrections. 

We can also investigate the q 2 dependence of the form factors fo* (? 2 )- It is known that the OPE for the correlation 
function is valid only at small momentum transfer region < q 2 < (mj, — m c ) 2 — 2Aqcd(w& — m c ). At the large 
momentum transfer region, we need to parameterize them in terms of phenomenological models. To achieve this goal 
we firstly analyze the form factors within the HQET framework, which works well for the b — > c transition. The 
matrix elements responsible for B — > Dq transition can be parameterized as |32j 

(D* Q + (P)\c W b\B(P + q)) = -i^m B m D , [ v +. (w)(v + v'), + Vd , (w)(v - «')„], (18) 
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TABLE I: Numbers of ff(0) and vt( w ) determined from the LCSR approach, where the uncertainties from the Borel mass, 
threshold value, quark masses and decay constants are combined together. For comparison, the results estimated in the QCDSR 
are also collected here. 
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where v — (P + q)/rriB and v' — P/mu* are the four- velocity vectors of B and Dg mesons, and w = v • v' = 
(m 2 B ± m 2 D , — q 2 )/2mBmD*- Combining Eqs. ((J) and (fT5)) . we have 

f?{q 2 ) = - 1 [{m Bi + m Di )r)f(w) - (m Bi - m Di )r)7 (w)], 
^/m Bi m Dz 

fni 2 ) = \l^[vtM + v7Ml (is) 

V m Bi 

with i = 1, 2 denotes strange and strangeless charmed scalar meson respectively. Similarly to the Isgur-Wise function 
£(v ■ v 1 ) for the s-wave transitions, heavy quark symmetry allows to relate the form factors rjf(w) and r/~ (w) to a 
universal function t 1 /2(w)[^] 

vtM + = -2t 1/2 (w), V f(w) - V -(w) = 2r 1/2 (w). (20) 

Different from the Isgur-Wise function t;{w), one can not employ the heavy quark symmetry to predict the normal- 
ization of Tj /ofi/;) [33j. 

Phenomenologically, one can parameterize the B — > Dq form factors in the small recoil region as 

r,t(w) = V ±(l) + af(w - 1) + bf (w - l) 2 , (21) 

The parameters rjf(l), af and bf can be determined by connecting the form factors derived in the LCSR and HQET 
approaches in the vicinity of region with q 2 ~ (mj — m c ) 2 — 2Aqcd(™& — fn c ). In this way, we can derive the results 
of form factors in the whole kinematical region, in Fig. (j4]) we take {q 2 ) as an example. The parameters related 
to all the form factors are tabulated in Table HI 

As discussed before, the power-suppressed form factors /~ in Table U suffer from sizable power corrections, which 
can even change the sign. Generally speaking, the corrections can be picked up by perform the heavy quark expansion 
of the current 

cT t b = c v ,Tib Vl - — !— c V2 Tii p± 2 b Vl + 7r~c V2 Tii p±ib Vl ± ... (22) 
Zm c Zrrib 

The last two terms in the above equation might give important contribution for finite quark mass, which could help 
to reduce the discrepancy among different approaches. In addition, the radiative correction may also help. 



IV. SEMILEPTONIC DECAYS 



The semileptonic decays B( s ) — > D^^lv are important measurements in the B factory which can be connected with 
the form factors directly. The differential decay width is given by: 
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(23) 
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FIG. 5: The q 2 dependence of differential decay width ^T(B° 
I — r (right figure). 



D*gl v{) for the final states with I = e, fx (left figure) and 



The q 2 dependence of these B® — > D*^l~Di partial decay rates are plotted in Fig. ([5]). Similar figures can also 
describe the B~ —J- Dq°1~Pi decays. The curve of the r final state is different from the light quark case for its mass 
effect. Integrating Eq. (j2"3"j) . we get the branching fractions of B^ — > D*^lu as grouped in Table UU The results 
from the constituent quark model, the QCD sum rules and the light quark LCSR are also listed here. Our result is 
slightly larger than the light quark LCSR as we have obtained large form factors. Note that the theoretical error is 
very large, which makes all the results are actually consistent. Besides, we can also find that the decay rates for the 
final state with r lepton are generally 3 — _4 times smaller than those for the muon case due to the suppression of 
phase spaces. The branching fractions for B° Dq° (2A00)lDi are also available, which is the first prediction for these 
decays, and we hope the future experiments can check our results. 



V. DISCUSSION AND CONCLUSION 



The charmed scalar meson spectroscopy has received many research interests both experimentally and theoretically. 
The -B( s ) decays provide ideal places to study their property. In this article, we employ the i?-meson light-cone sum 
rules to compute the B® — > L>*^~(2317) and B~ — > Dq°(2400) transition form factors at large recoil region, assuming 
L>*(f (2317) and L>o°(2400) being scalar quark-anti-quark states. With the help of HQET, we extrapolate the result 
to the whole momentum region, the q 2 dependence has beep plotted in FigQ. Our results are compared with the 
studies using the other nonperturbative methods, such as the light-quark LCSR, the QCD sum rules and the quark 
models. Considering large uncertainties, our results are consistent with these studies. Meanwhile, we also found that 
the power corrections should be large, which even change the sign of power-suppressed form factor /~ + . 

^so (2317) 

Subsequently, we utilize the form factors obtained using B-meson LCSR to estimate the semileptonic decays B® — > 
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TABLE II: Branching ratios for the semileptonic decays B® — ¥ Df (2317)Wi and B~ — > Dq°(2400)Wi with the form factors 
estimated in B-meson LCSR, where the results calculated in the conventional light meson LCSR, the constituent quark model 
and QCDSR are also displayed for comparison. 



B s D^ l vi I = e, I — t 

this work (6.0 ± 1.9) x 1CT 3 (8.2±^) x 1(T 4 

Light meson LCSR (2.3^;g) x 1(T 3 (5.7±|;|) x 10~ 4 

QCDSR[8] ~ 10~ 3 ~ 10~ 4 

Constituent Quark Model[6] (4.90 - 5.71) x 10 -3 

QCDSR in HQET[JJ (0.9 - 2.0) x 10~ 3 

B° — > Dq°1~ui I = e, [i I — r 
this work (8.7±|;g) x 10~ 3 (l.l±g;g) x 10~ 3 



D* + (2317)lDi and B~ — > Dq° (2400)iP;. It has been shown in this work that the branching fraction of the semileptonic 
—> D* + (2317)ZPj decay is around 6 x 10~ 3 for light leptons and 0.8 x 10 -3 for tau final state. The difference is 
due to the phase space suppression. The predicted values can confront with the future LHCb measurements. The 
predicted branching ration of B~ — y Z?q°(2400)/P; is slightly larger than — > D* + (2317)ZP;, and this observation 
can be tested at both LHCb and super B factories. 
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Appendix 



In the following we show the form factors from the 3-point B meson DA. 



f 3p 



f B (m c -m q ) [ p° mB ^ f°° de e -(^-« D j a )/M| 



^Id'TUd* Jo Jri mB-u £ 



xfM(m B Af + m c A± + Af - A± + Af ^+^7+"*^ ) 



V JO J V m B -Lj £ M B 



x(m B B 1 +m c B 2 + B 3 - B 4 + B 5 ^ T 



D 



(s -m D , 2 )/M% r i) m E 

doj I ^(Cf + Ct + Cf + Cf)}, (24) 



fir, \c m D* 2 >/ nl B pnamB roo jt 



2m B J Jrioms-u £ 



9 



where the functions Af(i = 1, 2, 8),Af(i = 1,2,..., 8) and Af(i = 1,2,..., 8) entering the integration are given 
below: 

A + , 2ao(2-37 lQ + (s a -q 2 )/m 2 B ) 1 + 2(s Q - q 2 )/m 2 B + 3(r c - ? 7o ) 

^1 = rj rj \{ipA-iPv) 

mBTJO TTlB 7/0 

, 6a »?o - 6r c 
= ml* 2 ^ 

m B rio 

A + = M-m 2 B f]$ + m 2 c + q 2 ) 
4 m B fjo 

At = *±?£-X A 

_ 2[m%(fj -2r c ) + ^-}(fjo+r c ) 

A 6 — To Xa 

m B r)o 

At = ~ 24a °f x A 
rriBVo 

,+ _ (-24m c )[m|(77 - 2r c ) + ^-^}(m + r c ) ^_ 

^8 — =a *A 

rriBVo 

r 2a (l-3r; + (so- g 2 )/m|) 2 + 2(a - g 2 )/m| + 3(r c - 7 ?0 ) 

A l = I — rj \WA-ipv) 

7775770 mBVo 

-6a T] Q - 6r c 

A 2 = j-r^ ipv 

A- = %*<l + *)x A 
m B rio 

a^mlvofjo + ^ml-^q 2 ) 

A 4 — To Xa 

m B % 

A 5 = — — X A 

mBVo 

_ 2[m%(r 10 -2r c ) + ^^](- Vo + r c ) 

^6 — =3 Xa 

m B r)o 

A f = ^^Xa 
._ _ (-24m c )[m|(7 ? o-2r c ) + ^-£](-7 ?0 + r c ) 

— Tg iA (25) 

777B 7?q 



B+ = 2a(2 - 37/ + (s - q 2 )/m 2 B ) - (1 + 2(s - q 2 )/m% + 3(r c - r/))(Vu - Vv) 

= (6a?7 — 6r c )ipv 

B+ = {-2am B )X A 

B+ = 2am B (-m 2 B f] 2 + m\ + q 2 )X A 

B+ = (f] + 2r c )m B )XA 

B+ = 2[m%(v - 2r c ) + g° 7 g2 ](77 + r c )X A 

B+ = (-24am B m 2 c )X A 

2 2 

B+ = (-24m e ) [77773(77 - 2r e ) + c _ q ](t? + r c )F A (26) 
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flf = 2a(l - 3?/ + (s - g 2 )/m|) - (2 + 2(s - g 2 )/™^ + 3(r c - r?))(Vu - W) 
B2 = (—Gar) — 6r c )'0y 
20(1 + 7?) ^ 

^3 = =9 — A A 

mBV 

D— o ^ 2 — . 1 + 7 7 2 V 2\\ v 

B 4 = 2am B (m B r)r) -\ — m — -q ))X A 

77 77 

S 5 " = (1 + V - 2r c )m B )X A 

2 2 

B 6 = 2[m%(fj - 2r c ) + ™ c ~ 9 ](-r? + r e )X A 
Sf = (-2iam B ml)X A 

2 2 

Sg = (-24m c ) [777^(77 - 2r c ) + TOc 7 9 ](r? + r e )y A (27) 



ar\ 77 



^ = ^[ 2 [m|(,-2 rc) + ^] ( , + , c) ^] 



C+ = ^[(-24^ 

C 8 + = A[(_24mJ[m|(77 - 27g + ^i^](77 + r c )^+ )y ' 4 
a?7 —77 ?; 

C A = —[20(771^7777 + -^m c - -g ) =3 — L, =7J0 

(177 77 77 77° 

6 = ^[ 2 [ TO sW-2r c ) + }(-V + r c )—r- 3 \ v = Vo 



3 Ji)=»7o 



C7 



dri [[ c) 773 lr ~ 



C 8 " = ^[(-24m c )[m|(r7-2r c ) + ^-^](-r7 + r c )^l^] ?)=7?0 (28) 

Where the notations 77 = cj + fa, /(r?) = (l + , X A (w,0 - / drX A (T,0, Y A (r,,0 = J drY A (r,t), and 

770 satisfies the equation 77S0 — (7777 + r^)m B + 77771^ = 0. 



B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 90, 242001 (2003) |arXiv:hep-ex/0304021] . 

K. Abe et al. [Belle Collaboration], Phys. Rev. D 69, 112002 (2004) [hep-ex/0307021|. 

J. M. Link et al. [FOCUS Collaboration], Phys. Lett. B 586, 11 (2004) [hep-ex/0312060|. 

B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 79, 112004 (2009) [arXiv:0901. 12911 [hep-ex]]. 

N. Isgur and M. B. Wise, Phys. Lett. B 232, 113 (1989); Phys. Lett. B 237, 527 (1990) ; Phys. Rev. D 43, 819 (1991). 

S. M. Zhao, X. Liu and S. J. Li, Eur. Phys. J. C 51, 601 (2007) |arXiv:hep-ph/0612008] . 

M. Q. Huang, Phys. Rev. D 69, 114015 (2004) [arXiv:hep-ph70 404032" . 

T. M. Aliev and M. Savci, Phys. Rev. D 73, 114010 (200 6) [arXiv:hep-ph/0604002| . 

K. Azizi, Nucl. Phys. B 801, 70 (2008) |arXiv:0805. 28021 [hep-ph] j. 

R. H. Li, C. D. Lu and H. Zou, Phys. Rev. D 78, 014018 (2008) [arXiv:0803. 10731 [hep-ph]]. 

S. Has himoto, A. X. El-Khadr a, A. S. Kronfeld, P. B. Mackenzie, S. M. Ryan and J. N. Simone, Phys. Rev. D 61, 014502 
(1999) [arXiv:hep-ph/9906376l . 

G. M. de Divitiis, E. Molinaro, R. Petronzio and N. Tantalo, Phys. Lett. B 655, 45 (2007) [arXiv:0707.0582l [hep-lat]]. 
G. M. de Divitiis, R. Petronzio and N. Tantalo, JHEP 0710, 062 (2007) [arX"iv:0707.0587l [hep-lat]]. 
R. -H. Li, C. -D. Lu and Y. -M. Wang, Phys. Rev. D 80, 014005 (2009) [arXiv :0905.3259 [hep-ph]]. 
I. I. Balitsky, V. M. Braun and A. V. Kolesnichenko, Nucl. Phys. B 312 (1989) 509. 
V. M. Braun and I. E. Filyanov, Z. Phys. C 44 (1989) 157. 



11 



V. L. Chernyak and I. R. Zhitnitsky, Nucl. Phys. B 345 (1990) 137. 

M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147 (1979) 385; Nucl. Phys. B 147 (1979) 519; Nucl. 
Phys. B 147 (1979) 448. 



A. Khodjamirian T. Mannel and N. Offen, Phys. Lett. B 620 (2005) 52 [hep-ph/0504091| ; Phys. Rev. D 75 (2007) 054013 
[hep-ph/06lTl93] . 

A. G. Grozin and M. Neubert, Phys. Rev. D 55 (1997) 272. 

A. G. Grozin, Int. J. Mod. Phys. A 20, 7451 (2005). ' 

Y. -M. Wang, Y. -L. Shen and C. -D. Lu, Phys. Rev. D 80 (2009) 074012 [arXiv:0907.4008l [hep-ph]] . 

A. Khodjamirian, T. .Mannel, A. A. Pivovarov and Y. -M. Wang, JHEP 1009 (2010) 089 [arXiv: 1006.49451 [hep-ph]]. 

I.I. Balitsky and V.M. Braun, Nucl. Phys. B311,541(1989). 

J. Beringer et al. (Particle Data Group), Phys. Rev. D86, 010001 (2012) 

P. Colangelo, F. De Fazio and A. Ozpineci, Phys. Rev. D 72, 074004 (2005) [arXiv:hep-ph/0505195| . 
T. M. Aliev and V. L. Eletsky, Sov. J. Nucl. Phys. 38, 936 (1983) [Yad. Fiz. 38, 1537 (1983)]. 

P. Colangelo and A. Khodjamirian, In Shifman, M. (ed.): At the frontier of particle physics, vol. 3, 1495-1576 
|hep-ph/0010T75] . 

H.G. Dosch, E.M. Ferreira, F.S. Navarra and M. Nielsen, Phys. Rev. D 65 (2002) 114002 [arXiv:hep-ph/0203225j . 



Y. M. Wang, H. Zou, Z. T. Wei, X. Q. Li and C. D. Lu, Eur. Phys. J. C 54 (2008) 107 |arXiv:0707.1138l [hep-ph]]; Eur 
Phys. J. C 55 (2008) 607 |arXiv:0802. 27621 [hep-ph]] . 
P. Colangelo, G. Nardulli and N. Paver, Z. Phys. C 57, 43 (1993). 

T. Kurimoto, H. n. Li and A. I. Sanda, Phys. Rev. D 67, 054028 (2003) [arXiv:hep-ph/0210289 . 



F. De Fazio, In Shifman, M. (ed.): At the frontier of particle physics, vol. 3, 1671-1717 |hep-ph/0010007 



